Structures of quantum 2D electron— hole plasmas 



V S Filinov\ M Bonitz^, H Fehske^, P R Levashov\ V E 
Fortov^ 

^ Joint Institute for High Temperatures, Russian Academy of Sciences, 
Izhorskaya 13 bldg 2, Moscow 125412, Russia 

^ Christian- Albrechts-Universitat zu Kiel, Institut fiir Theoretische Physik und 
Astrophysik, Leibnizstrasse 15, 24098 Kiel, Germany 
•^Institut fiir Physik, Ernst-Moritz-Arndt-Universitat Greifswald, 
Felix-Hausdorff-Str. 6, D-17489 Greifswald, Germany 

Abstract. We investigate structures of 2D quantum electron-hole (e-h) plasmas 
by the direct path integral Monte Carlo method (PIMC) in a wide range of 
temperature, density and hole-to-electron mass ratio. Our simulation includes a 
region of appearance and decay of the bound states (excitons and biexcitons) , the 
Mott transition from the neutral e-h plasma to metallic-like clusters, formation 
from clusters the hexatic-like liquid and formation of the crystal-like lattice. 



PACS numbers: 71. 23. An, 71.55.Jv, 52.65.Pp 

1. Introduction 

Strongly correlated two-dimensional quantum Coulomb systems are the subject of 
intensive discussions. In particular, it is known that the competition between 
electrostatic and kinetic energy in an electron gas may be the reason of an unusual 
phase diagram of a 2D system of electrons. The liquid state of such system is 
stable when the kinetic energy dominates while the electrostatically favored "Wigner" 
triangular crystal is stable in the opposite case. If there is a strong competition 
between these two kinds of energy, different situations are possible. The question 
under discussion is the existence of the intermediate anisotropic liquid phase (hexatic) 
under melting of crystal into isotropic liquid. Moreover, the physical mechanism of 
melting can be influenced by the interaction with substrate and defects. Currently all 
mentioned phenomena can be extensively investigated by a consistent first-principle 
numerical simulation, and in this brief paper we present the most interesting results 
of our numerical experiments by the direct PIMC method. 

2. Path integral Monte Carlo approach 

Let us consider a two-component neutral e-h plasma. Thermodynamic properties 
of such plasma are defined by the partition function Z, which for the case of 
electrons and Nh holes {N,, = Nh), is given by Z{Ne, N,,, V, /3) = Q{Ne, Nt, f3)/N^\Nh, 
with Q{Ne,Nh,l3) = /y ci'^ c;<7 ^, cr; where /3 = l/ksT, a = {ae,<Jh), r, a^. 
denotes the space and spin electron coordinates, while q, denotes the space and 
spin hole coordinates. These expressions are valid for ID, 2D and 3D cases. The 
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pair distribution functions for a binary mixture of quantum electrons and holes can 
be written in the form gab{Ri - R2) = gab{Ri,R2) = Eo- !v ^1^'^ ^(^^ ^ '^(-^2 - 
Q2)p{q,r,a,(3)/Q{Ne,Nh,f3), where a and b label the particle species, i.e. a,b — e 
(electrons), h (holes) and Q,R denote the two-dimensional vectors of coordinates, 
Qi 2 ~ ^1,2 E^nd 2 — 91,2- The exact density matrix p{q, r, a, f3) of a quantum system 
for low temperature and high density is in general not known but can be constructed 
using a path integral representation [1,3]. Thus we take into account interaction, 
exchange (through permutation operators) and spin effects both for electrons and 
holes. This procedure gives the expression for the density matrix as a multiple integral 
which is well suited for efhcient numerical evaluation by Monte Carlo techniques, 
e.g. [1]. We mention here that for all obtained results, the maximum statistical and 
systematic errors is not more than 5%. 

3. Simulation results 

We are interested in strong Coulomb correlation effects such as bound state (excitons, 
bi-excitons, many particle clusters), their transformation by the surrounding plasma 
and their eventual breakup at high densities (Mott effect). Beyond the Mott density, 
we expect the possibility of the hole crystallization if the hole mass is sufficiently 
large [2]. Below, the density of the two-component plasma is characterized by the 
Brueckner parameter rg defined as the ratio of the mean distance between particles 
d — [l/irine + 7T./i)]^^^ and the 3D exciton Bohr radius as, where and Uh are the 
electron and hole 2D densities. The dimensionless temperature will be presented as 
a ratio of the temperature and the 3D electron- hole binding energy (Rydberg) , which 
includes the reduced effective mass and dielectric constant. 

We analyze some spacial distribution functions and related spin-resolved typical 
"snapshots" of the e-h state in the simulation box for different particle densities, 
temperatures and hole-to-electron mass ratios M . According to the path integral 
representation of the density matrix, each electron and hole is represented by several 
tens of points ( "beads" ) . The spatial distribution of the beads of each quantum particle 
is proportional to its spatial probability distribution. Figure [T] shows that, the typical 
size of the cloud of beads for electrons is several times larger than the one for the 
heavy holes. At low temperature and low (ts — 6) and middle (r^ — 2) densities 
practically all holes are closely covered by electron beads. From a physical point of 
view this means that electrons and holes form bound states, i.e. excitons, bi-excitons 
and many particle clusters. The existence of the bound states is also supported by 
the behavior of the pair distribution functions, exhibiting pronounced maxima at the 
distances of about half and a bit more of the Bohr radius. We note that raising the 
temperature at a fixed density leads to a temperature-induced ionization of the bound 
states. As a result we found a substantial number of free electrons and holes. 

From figure [1] it follows that the growth of density results in the increase of the 
number of particles in clusters. The structural analysis of large many-particle clusters 
show the hexagonal ordering of heavy holes inside liquid-like clusters (lower right 
panel). In this panel there are only two clusters: one is in the center of the Monte 
Carlo cell, while the second one is divided in four parts (in each corner) due to the 
periodic boundary conditions of the Monte Carlo cell. Here besides the inner normal 
hexagonal structure the holes due to the strong Coulomb repulsion form the filament- 
like structures of the clouds of beads at the bounds of clusters. The filament-like 
structures are fully 2D topological effect as in 3D case the holes in analogous many 




Figure 1. (Color online) Pair distribution functions (left column) and snapshots 
of the Monte Carlo cell (right column). Left column: g^e (black solid line), 
ghh (red dashed line), jeh (blue dot-dashed line) at temperature T/Ry = 0.007, 
the hole to electron mass ratio M = 800 and densities related to the Brueckner 
parameter = 6 and = 2. Right column: red and magenta clouds — holes, 
yellow and cyan clouds — electrons with the opposite spin directions. 



particles clusters have normal liquid- like ordering [3,4]. 

If the particle density is high enough (r^ — 0.25) the electron wavelength becomes 
larger than the mean inter-particle distance d and even larger than the size of the 
Monte Carlo cell which is seen by the large extension of the clouds of electron beads. 
For < 0.5 clusters become unstable because two electrons bounded to neighboring 
holes start to overlap allowing for electron tunneling from one cluster to the other 
(Mott effect). Since the hole wave-length is significantly smaller than the electron 
wavelength, it may still be smaller than d and the structure of the hole beads resembles 
a liquid-like or crystal-like state. If the hole mass exceeds a critical value, the holes 
may even form a crystal- like structure [2] (lower right panel in figure [2]). Here the 
holes form a crystal-like structure, while the electron density demonstrates the Bloch 
oscillations. At this very high density the type of the hole 'crystal' is influenced by the 
boundary conditions of the Monte Carlo cell (finite size effect). The detailed analysis 
of the type of the crystal-like structure is possible for a much more bigger number of 
particles in the Monte Carlo cell. 

4. Conclusion 

In this paper we have presented a computer simulation analysis of strong Coulomb 
correlations in dense two-dimensional two-component quantum plasmas at low 
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Figure 2. (Color online) Pair distribution functions (left column) and snapshots 
of the Monte Carlo cell (right column) at temperature T/Ry = 0.007, the hole- 
to-electron mass ratio M = 800 and densities related to the Brueckner parameter 
= 0.5 and = 0.25. Notations are the same as in figure [T] 



temperatures. In particular, the formation and dissociation of bound states, such 
as excitons, bi-excitons and many particle clusters, is analyzed and the density- 
temperature regions of their occurrence is identified. We have shown that, above 
the Mott point, two-component plasmas with large mass anisotropy show interesting 
Coulomb correlation phenomena: with increasing density holes can undergo a phase 
transition to a Coulomb hexatic-like liquid and to a Wigner crystal which are 
embedded into a degenerate electron gas. The crystal-like formation in a two- 
component plasma is possible for large enough hole-to-electron mass ratios [2]. 
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